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ABSTRACT 


A computer program for calculating internal supersonic flow fields with 
chemical reactions and shock waves typical of supersonic combustion chambers 
with either wall or mid-stream injectors is described. The usefulness and 
limitations of the program are indicated. The program manual and listing are 
presented along with a sample calculation. 


The work reported herein was supported by the National Aeronautics and 
Space Administration under Grant NGR-33-016-131 . 

t Assistant Research Scientist, New York University Aoiospao.' Labor: tory 

+f Research Scientist, New York University Aerospace Laboratory 


NASA-CE- 112303) DESCRIPTION OF A 
COMPUTER PROGRAM TO CALCULATE REACTING 
S U PE RSO NIC' 'INTER N A L FLOW FIELDS WITH 
SHOCK HAVES USING VISCOUS (New York 
Univ.) 100 p HC $7.00 CSCL 20D 


N73-?0311 


Unclas 

G3/ 1 2 67736 



TABLE OF CONTENTS 


ABSTRACT 

INTRODUCTION 1 

ANALYSIS 2 

DESCRIPTION OF COMPUTER PROGRAM 3 

PROGRAM MANUAL 4 

SAMPLE CASE AND RESULTS 5 

CONCLUSIONS 5 

REFERENCES 7 

TABLE 1 
TABLE 2 


FIGURES 



LIST OF FIGURES 


FIGURE 


la,b 

2 


3a 

3b 

4 

5a-d 


6a -i 


7 


Simplified scramjet fuel injector arrangement 

Simplified combustor entrance flow field regions 
and nomenclature 

Fuel injector geometry 

Local injector wall angle 

Initial flow profiles 

Flow properties along injector wall 

Flow profiles at various axial stations 

Characteristic mesh 


PAGE 


72 

73 

74 

75 

76 

78-81 

82-99 

100 


1 

2 

3 


TABLES 

Program listing 
Input cards 
Sample input data 


8 

62 

69 



INTRODUCTION 


The structure of the flow field within a supersonic combustion chamber 
designed to operate over a range of flight Mach numbers is complex due to many 
interacting phenomena. This flow field is permeated both by nonuniformities 
and discontinuities. The nonuniformities inherent in the incoming stream are 
produced by the inlet flow, fuel injectors, and combustor walls. The non- 
uniformities exhibit themselves in a variety of ways among which are: 

1) Nonuniform initial conditions such as flow direction, Mach 
number, and total pressure distributions. 

2) The boundary layer on the inlet walls introduce shear layers 
with very large gradients normal to the streamlines and 
possible reverse flow if separation occurs. 

3) Oblique shock intersections produce slipstreams as well as 
rotational flow and vortical layers. 

4) The fuel injector and combustion phenomena produce non- 
uniformities in the gas composition and phase. 

5) Other nonuniformities are due to three-dimensional fields. 

In view of the above complexities, a truly realistic theoretical descrip- 
tion of a supersonic combustion flow field presents a formidable task. A 
first attack at this task is attempted here for a more simiplified flow field. 
The assumptions inherent in the present analysis are a steady, viscous, two- 
dimensional (or axisymmetric) , completely supersonic flow field. In view of 
the last assumption, the boundary layer formed near solid walls is disregarded 
here. Also, changes in the local sound speed produced by either variation of 
species, or temperature rise as produced in exothermic reactions, or excessive. 



flow deceleration produced by pressure rises so as to drive the flow subsonic 
are not allowed. The steadiness of the flow is imposed on the fluid mechanical 
time scale as well as on the chemical time scale. Therefore, combustion 
unsteadiness, "sputtering" as produced by conditions in near flamability limits 


nharmodynauric conditions of the fluid 


is not considered hare. 


The chemistry 


is assumed to occur at a finite rate or frozen and the reactions are assumed 
to occur between air and gaseous hydrogen. The viscous transport phenomena 
can either be laminar or turbulent. 

The immediate problem at hand is to predict the flow downstream of a 
given initial station and bounded by an upper and lower wall with two either 
converging or diverging shocks and with chemical reactions as shorn on Figs, la 
and lb. Region one can be either uniform or nonuniform. The flow properties 
across the shocks are calculated assuming no gas composition change across 
them, i.e. as in an ideal gas using the local value of the ratio of specific 
heats . 


ANALYSIS 

Since the flow field is assumed to be completely supersonic (both air 
stream and fuel stream) the field can be considered to be "quasi-hyperbolic" 
if viscous stress terms due to gradients along the streamlines are neglected 
and the viscous stress terms due to gradients normal to streamlines are 
approximated by constants in a small neighborhood. Under these assumptions 
the governing equations can be cast in characteristic form with driving-functions 
on the right hand side. The driving functions depend solely on the viscous 
transport and chemical phenomena. A derivation of the compatibility equations 
along the characteristic lines may be found in Refs. 1, 2, and 3 along with the 
chemical reactions, reaction rates, and calculation procedure. The imposition 
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of the proper boundary condition is observed in the present computer program. 
The turbulent viscosity model used in the present program is 

U . = K r, o u 
M t % • o o 


where r . is a measure of the size of the eddies and is determined as the 

'2 


height where D u = • % (n u +r. u ) 

^ v e e ^o o 

o denotes centerline or wall conditions. 


where e denotes external conditions and 


COMPUTER PROGRAM 

The computer program proceeds from a given initial data line at a constant 
axial station to the next station by taking a calculated step size satisfying 
both parabolic and character istic as well as chemical stability criteria. A 
listing of the program is given in Table 1. 

In the main part of the program the initial profile is set up and the type 
of, flow (i.e.j converging or diverging shocks) is selected. The main loop for 
each step begins at statement number 6789 and ends at the statement prior to 
1572. The following is a list and description of the sub-routines in the pro- 


gram: 

MESH 

SHKINI 

HERMAN 

CLEM 

SOLT 

ETHANE (T) 
PROP(T) 
ETH 2 (T) 


traces a characteristic line 

determines the intersection point of converging shocks 
solves for chemical time 


fudges the chemistry for ethane-air reaction 
fudges the chemistry for propane-air reaction 





SURFAC ( IND ,L , X , R , DR) 
CUBIC (C,Z) 

COEF (I ,T) 

THERM ( IND ,L) 


reads in wall slopes and computes wall position 
solves for shock wave angle 

gives thermodynamic coefficients, constant, fits 
calculates enthalpy 


UFSC (THS) 
STEP(DELX,N1,NP2) 

INI ER ( RAT , KL , N1 , N2 ) 
FIND 

ERROR(IIII) 

SHOCK 

FLIP(IMB ,N1 ,NPTS ,XMSR) 
FLOW 

POINT ( ID , NPO , IU , IFM) 
BODY ( JXT) 

FOCUS (TI , PRESSI , RKOI . ) 
CHEMP ( IWD , DX , L ) 

COMPS (IND) 


calculates upstream £iov? if nonuniforrn 
calculates step size 
interpolation subroutine 

locates intersection point of streamline with given 
data line 

writes out error message number 
calculates shock point 

stores new station into old - updates data 

calculates an intergral point 

adds internal point near shock 

calculates wall point 

calculates finite rate chemistry 

calculates chemical production term 

calculates viscous dissipation terms (viscosity) 


PROGRAM MANUAL 

The program can be used to calculate the two flow fields described in 
Figs. 2a and 2b. The control cards and input cards necessary are given in Table 
2. A sample input data is shown in Table 3. The present program is not 
equipped to handle a slip stream, discontinuity as would occur immediately at the 
lip of the injector. The initial profile data must be continuous. The slip- 
stream can therefore be approximated by a strong gradient in the y direction. 
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Sample Calculation - Hall Injector 


A sample calculation using this program was made for the wall injector shown 
in Fig. 3. The H ? jet was chosen 0.02 ft. high. The jet initial Mach number is 
2.0 and the air external stream is 2.7. Other initial conditions (temperature, 
velocity, pressure, Mach number, flow direction, and ^ mass fraction) are shown 
in Fig. 4. The fuel injector shock, was located away from the combustion zone in 
this calculation. 

The lower wall geometry has a 10° slope discontinuity (expansion) at the 
mouth of the injector. The expansion wave from this point propagates into the 
combustion zone. The wall angle distribution for this injector is shown in 
Fig. 3b. The wall is assumed inviscid (i.e. no boundary layer) in this calcula- 
tion. The pressure, Mach number, velocity, and temperature distributions along 
the wall of the combustor are shown in Figs. 5a through 5d . The corresponding 
distributions calculated assuming Prandtl-Meyer expansion of the flow along this 
wall are shown for comparison. The differences are due to down running waves 
from the combustion zone and the external stream. Flow profiles at downstream 
stations are shown in Fig. 6a through 6i. The temperature profiles show a 
characteristic peek typical of a flame. While the pressure profiles show the 
propagation of combustion induced waves . The species profile show the diffusion 
of the fuel into the flame water from the flame. 

A plot of the characteristic and streamlines through the combustion zone is 
shown in Fig. 7. The flame zone. edge of mixing and other flow characteristics ar 
shown for comparison. 


RESULTS 

A useful tool for analyzing a simplified supersonic flow with chemical 
reactions is presented. In spite of its Limitations, the computer program can 
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can be used 1) in assisting in the design of fuel injectors and supersonic com- 
bustors, 2) in segregating the highly coupled fluid -mechanical effects of com- 
bustion and mixing in supersonic streams, and 3) in evaluating turbulent eddy 
viscosity models by comparing experiment with theory. 
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Table 1 Program Listing 


rOi'-' a USTOR 


J 

J 

U 

O 

o 

o 

o 

o 

u 

o 


o 


o 
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N 

‘O. 
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0 i b'FK S 1 ON 5,v A ! 3 ( 1 0 ) « Y ( y ) , Y V ( 7 ) 

T j v - = r> t 
T .V X = T T M - -if * ;> 

YN f \ ) =Y ( ] ) -f SiMALr- ( ] )*7 ; M^Syi.LP ( 8 ) *TNX 
Yv ( ? ) =Y ( 2 ) +SVAL5 ( 2 ) -* T ! *•* F + S A L = . ( 6 ) * T X 

Y N ( A ) =Y ( 6 ) -Sf-'ALB ( ? ) *1 ? y “ + Sf-'A LR ( ) * T,\X 

YN ( ”=> ) =Y ( -a ) -i- S M A L B ( u )-ifT IMF + 5MALB ( 3 ) *T\iX 
YM ( A ) = CC — ( Ym ( 1 ) + YiY ( a ) 5 / ? . 0-YN ( “ ) 

YM < 8 ) = R P - ( Y\i { 2 ) + V K ( R ) i-Ylx (?) ) / p „ 0 

return 

p,\!R 


SURROLJT T V; P - IMAK'- ( T ). 

T = ], . / ( 1 . 108 7/T-.09A97 ) 
RETURN 
p Ni r> 


SI )a !? o i j T t m p ppno( i i 
T = 1 , / ( . 'V86/T + . 238 1 ) 
RETURN 

tr Nj on 


c;. p oOl I T t V> c r T u ? ( Tj 
T=1 . 241+.05524*T 

return 
p M °i 
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-o 

o 

o 


Sji5 POO [ T \f . g .Un)--Ar ( I NO , I_ » X , R OR ) 

D I MENS] ON XRC( 100 > ,RRO( ]C0J »TMC( • 00) » THA { 100 ) ,?7 { i 00 > ,P { 1 CO ) 

rv/o.017^5 

T=L _ 

r , 0 TO ( 6 , 7 ) » ! NO 
-> T = L-2 

6 r-0 T 0 ( ] > ? » ^ > A i t j 
]. R- AO ( F • i 3? ) C 3 . i COWL 
J =1 
jj=rp 

n O K=liS 0 

D ( '< ) = 0 . 

1 R PT (■< ) = 0 . 

.R E A 0 (C..10?) (XRCJO »< = J»JJ) 

RFAO (F,}02) PPC(Ji 
5 O 0 Q <= J , J J 
Tua ("< )=Tur(<) 

Twr ( .< ) = T AM ( T Hf ( < l > : rv ) 

I F ( K o F a . ' . OR . K . c O . = 1 ! GO To Q 

RRC ( K. ) =0 . s# < THC « < ) +THC ( K-l i ) * ! XRC ( !< » -XRC ( .<- 1 ) ) +RRC ( :<— 1 ) 
q r O N T T N U r 

IF ( I . ~Q. 0) .GO TQ jo 

WRtTF (6,]03i 

Vi‘R tTF ( f > » 105 ) ( XRC ( J ) * RRr ( J ) ♦ THA ( J ) » J= i ♦ J J ) 

10?. FORMAT ( 55X»?=K CENTER BODY CO-ORDINATES // »4SX»2H X » 14X» 

]].?X) 6 K TH r TA / ) 

T 05 FORMAT ( 4 0 X * 5 ( ? X E 1 4 . 6 ) ) 

RETURN 

19 X = XRr ( 51 ) 

R = R P r ( =, i 5 
WRIT- ( 6 , 1 0 A ) 

C WRtTF ( 6,1 05) (XRC ( J) »;RRC (J) » D CJ ) * J= 5 1 > J J ) 

W R T i F ( 6 » 1 0 5 ) ( X R C ( J ) > R R C > J 5 > i HA' J) » J = 5 1 » J J ) 

1.04 FORMAT ( ]ril ,67X,24H COWL CO-ORDINATES //»48*»2H X»14X» 

1 1 OX t 8H T H F T A /) 

R c TURN 

2 J = ? 
r p = i 
JJ = CR 

13 00 1.0 < = J,JJ 


M = < 


]F (X „LF, XRC(K)) GO TO 12 
] 0 oon'T i f'.’ijc 

I ? 00= ( T HC ( M ) - 7 H.C ( N- ] ) ) / ( XRC ( M ) -XRC ( V— ' ) ) 
po 1 - X-XRC ( M-l ) 

R = . =. * D 0 * D D 1 * * 2 - D 0 1 T K C ( ■•' - 1 i + R R C ( M - 1 ) 


PR = DO-k-DI i + ThC ( M- i ) 

PS = p ( (■■'- l ) ~r>01. / ( XRr ( V ) -XRr ( M- 
R 0 = ° T ( M-i j+DDl / ( XRr ( y ) -XRC ( v- 

R r TIJRN 


J = s 

T 

JJ = 

n 0 # 0 4- 

R F a 0 

< 5 , 1 0 ? ) 

RF Af) 

( F , 1 0 2 ) 

REA D 

( 3 » ] 0 2 } 

° F A D 

( 5 , 1 0? ) 

R f a 0 

( F , i 02 ) 

r-0 T 0 



OWL 

( XRC U ) ».<=J> JJ) 
RRC ( J ) 

( iKC( | <)*'<=J»JJ) 
( pT(k)»k=j»jj: 
( p ( y. ) » k = j , J j ) 


} ) 


)-? ( 
) -pi ( 


1 ) ) 
l > ) 


r\ 


r 

r 
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/' 


O 


O 


O 


4 j = 

jj = ".o.o + rowL 

r-C) TO 1 1 

10? POOKA T ( 7C1 0-.6 ) 
c-^ri 


U 


O 


o 


O 


O 

o 

o 

o 


o 

o 

o 


SUFROUT t rp rUPTCtr-Z ) 
p j v, -MSi 0 N C ( 4 ) 

acos ( x ) = aTa.n ( sort ( i . o-x **2 > / x ) 

. D = -c { R >**?/?. 0 + c ( 7 ) 

0 = ? . 0 * f ( - } -:;-*•}/? 7 . C - C ( 2 ) *C ( ^ ) / 3 . 0 + CM) 

rso = - 0 . 5 * 0 / sort (-p*«--s/ 2 7 .o> 

IF (ABS(RSO) . lE. 1.0) GC TO 1 
7 = 0. 

R-TL'RN 

1 Pht =AfOS ( RSO) 

. - p ( P H T } o 0 0 * o 0 3. ■» 9 0 3. 

9 00 PHT = 3 . 141 593 + PhI 
oOl T"M=R .Oh;- SORT ( -P/9 .0 ) 

XI = T FM*COS ( PH T /3.0 1 

Tfr ; v = PHT/3.0 + ' 2 .09439510 

X?= T cw*fOS ( T “PM ) 

Tcrv. = PHT./3.0 + 4.18879020 

X 2 = T- v*rOS ( TRRV, } 

TP ( >; ? - X 3 ) 150*150*160 
150 Y 1 = & M A X 1 (XI * X 2 ) 

Y ] = A >i t m 1 I Y 1 * X 9 ) 

r-0 T 0 3 5’ 

160 Y 1 = A T N 1 l X 1 * X 2 ) 

Y 1 = A :X A X ] . ( Y 1 * X 3 ) 

1 75 Yi=Yl -r { 3 ) /3 .0 

8 .7 = Y 1 
9 p T U 9. M 

FMPl 


o 

0 


o 


G 


O 
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3 


o 


u 

SU" 

p 0 ■> j T T M p l T * T 


COv 

•V{l<\i/COr/ A ,p ,c«r>» 


I- ( 

7-1 ODD. ) 1.0 * 10 * 20 


10 r, 0 

TO i 1 1 1 1 & > P > i 1 < l 

• 

11 A 

= 7 , l 460 64 0 p- Oft 


F 

= 4 . 1 03 ? 1 lA^-O^ 

0 

r 

= -0 , 6 1 1 F 3 3 2 -- 06 

- 

n 

= G . 3 ] ?26f ?F-na 


p 

= -3.309 3421- - I? 


p 

=-?. 67253^24 02 


r 

= -1 pippoOF Co 


r.n 

T 0 40 

o 

12 A 

= 3.7159945F Oft 


C 

= -2 . 9 1 67?.=,ft“-0P 


r 

= 8. . 963 7 2 5 3 - — 0 4 

o 

n 

= -6.2 7 9 6 7 I 6 ~ - 09 


rr 

= 2 . 7 0 8 2 1 3 0 “ - I ? 


p' 

=-1.0576706- 05 

o 

r~ 

= 3 . 9 0 fi 0 7 0 4 p Op 


no 

T Li A 0 


13 A 

= 4. 1 5650] 6 r Op 

o 

D 

= -] „ 724445A.P-0^ 


r 

— p # 6 9 ft 2. ^ 1 fi p* ~ 0 A 


n 

=-4 . 5Q30044F-0° 

o 

r 

= 1 .4233694-- 12 


rj- 

= -?«0?8f?77OF 04 

o 

,c 

= — 6 , 8 6 1 6 2 4 6 p - 0 1 

r-n 

TO 4 0 

* 

14 ft 

= 3 . 6 9 1 6 1 4 P c OP 

0 


- - ]. . 3 3 3 2 5 6 2 F - 0 ^ 

r 

= 2,4503 j 0 0 “ - 0 6 

• 

ro 

=-9. 7688 341 p-1o 

r*\ 

c ■ 

-—9 . 9772234F— 14 


p 

= - 1 , 062 8 3 36 c 03 


rz 

= 2.2874980" 00 


r- 0 

TO 40 

o 

1 3 A 

= 2 . 5000000- ’00 


p 

= 0.0 


r 

= 0.0 

o 

n 

= 0.0 


r 

= 0 . 0 


c 

= 2 . 3 4 7 0 4 9 7 c 04 

o 

r. 

= — 4 . 600 1096- -01 


r- 0 

T 0 4'0 


16 ft 

= 3.0218894- 00 

o 

Q 

=-2 . 1 73724pp-Op 



= 3 . 754220PP-0F 


r> 

r-2.944 7200F-CP 

.0 

c 

- O.0777547P— 17 


P 

= 2.9137190- 04 

o 

r\ 

= 2.6460.076- Op 

r-.n 

TO 4 0 


17 ft 

= 7 .8234 7 OR c OP 


P 

= - 1 1 1 R 7 7 2 o: p - P 0 

o 

r 

'= 1.24663199-06 


n 

= - 0 IP''- O p O A p — 1 pi 


C* 

= -5.23 4 6 5 5 .1 F - 1 4 

o 

r 

= p 5 0 4 7 7 8 7 r n ? 


24 



o 


o 


n 

= 3.825^029^-01 



r.Q 

7 0 9 0 

o 

20 

c, 0 

TO (23,26^2°. »21»27 

21 

A 

= 2.C43S697P Co 



D 

= 9.1 1371 10- -04 

. 


r~ 

= - 7 . 3 9 c 3 = 5 1 p _ 0 c 

o 


n 

=-2 . O 33 1907- -1 1 



c 

= 2 . - 3 9 3 7 9 1 r — 1 6 




= - P. .940100 2 r Op 

o 


r_ 

= — 1 .6961 3 ? 9 c 0 ^ 



r-O 

TO 4 0 

o 

22 

A, 

= p. c 976129p Op 


O 

— 7 . 814-3 6 03 *- — 0 a 



r 

= - 2. 23 3 997 op- 0 7 

o 



= 4 . 2990 1 59- —1 i 


c 

= -3 . 3^60204p-1 =, 



c 

=-i.L?2 7 9i8P 0^>. 

o 


r. 

= Oo 


r.O 

T 0 4 0 


23 

A 

= 2. 6 70 7 5 9 2 r Oo 

o 


0 

= 7 . 0 3 1 7 1 1 9 p - 0 0 


r 

= -B .9361 5 7 0 p - 0 7 



r> 

= I . ] 7 9 0 8 9 3 9 — Vo 

o 


c 

= -6 . 1 9 7 3 5 6 8 9 - I 9 


r 

= -2.9 R P 8 9 9 4 p 09. 



r: 

= 6.8838391 F Oo 

0 


r -0 

TO 40 

24 

A 

= 2.8949761P 0 n 



P 

= 1 . 99 76 3 16P-03 

- 


r 

=- 6 . 2969?9^P-0 7 

o 


r> 

= 1 , J 3 1 3 R 4 9 P - 1 0 



c 

=-7.6R97C 7 0 c -l9 

* 


r 

= — R . 9 0 1 7 4 9 9 p j- 0 2 

o 


n 

= 9 . 0 0 0 2 8 7 ci p On 



r 0 

TO ^0 


23 

A 

= 2.30000009 CO 

o 


P 

= 0.0 



r 

. = 0.0 




= 0.0 

o 


c 

= 0.0 



p 

= 2.947049 7 - 04 

o 


r 

=-4.60010963-01 


r.O 

TO 4 0 
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A 

= 2. 93725 6 7c Oo 

o 


p 

= -i . 64221 9 0 p - 0 5 


r 

=-P. 601792 1F-0G 



n 

= 5 . 9643621 F-l 2 

o 


p 

= — “ . 9 7 4 3 6 0 8 p - i 6 


cr 

= 2. 9.2 2 0007 p 04 

• 


c 

= 4.9467942P Oo 



r -0 

TO 4 0 

.0 

2 7 

A 

= 2. 8 S 5 5 9 49 P Go 



' P 

= 9 . 983 5 06 i p-04 



r 

--2.2 879904^-07 

0 


. 0 

•= 1.98027653-11 



r 

= — 3 . < 4529403-18 



c: 

= 3. fi.« 1 1792" 03 

0 


r 

r '.. 5 59 7 0 ] Pr r,n 


4 0 

D : ' 

■| 1 

r~\ 


r m 

- 

'■-.J 






J 

3 

■“N 

K 

D 

*\ 


S~\ 

ys 


O 

O 

o 

o 


SUQ vpi.i i i NT * ' O S r { Tjj^j 


CG V M O N / u l .< 7 / X MU ( i 
- \ / 


■ xy.'jN 


: .0 ! 


C ) , E vC T ( 7 • 3 = U ) 

J , VSN ( ■».) ; YS ( 3 ) 


•/.i5u i 
T ;-:3 

COMHG::/EL<a/ PHCSj>!( ~ ) »~ V .SN( ■> ! ,ASn; ( > ~-5 « T S.X ( 7 ) <X 

C 0 V M G *•:/“. L K t A / X '•* 'J S i\ ( ' ] • XMU3 i * ) ,kTS\ i -s ) > S iG'MA < "C ) 

COMMON /=L <6 / Hi N ( I = n ) , { •; =0 ) »YP(o) »Y( 1 f= 0 5 » Y :v ( 1 

CO v,MO.*: / hL / TH < 1 = 0 ) ,7 HN ( 1 ' 0 } 5 G i 1 O') » G ( 1 = 3 ) » T i 1 


CO'-' 

\ \ ,' *• : . 
i i ■. 

/ 

- L < 1 

7 / H X ( 3 

COM 

.* .**' ; 

/ 

c j_ ] 

'■ / O ( 1 

COM 

^0\ ! 

/ 

E L -K 1 

r / R ( } 

COv 

v * *' j 

/ 

-■ .< 1 

n / F i ( 

cn \' 

/ p, f 2 

/ 

- 1_ :< 1 

Q / m ? D 

COM 

•'Gt', 

/ 

L < ? 

0/ N'LO 

COM 

* /“■. •. 

' ’ i S. 

/ 

El <7 

1 / JP, 

COM 

v n n s 

/ 

EL <7 

7 / r I N » 

COM 

MOM 

/ 

/ 

UP] / 

GAMAU 

COM 

MON 

/ 

U?? /XV: IC UP ( 

X M 1 

; ex 

o t 

i OXA 

,0X7] 

Xy-a 

( oz 

■* 

•* Ql7 

♦ 0 z 7 , o 


-iXW (I'C; , a L ? ( 7 


. M ) , A. L r" N ( v » 2. 


A ; 1 V { *• L i \ !■ - ^ *> .■ 1 j r^nU:\ i i A C ) 

- C) < 5 7 -\- [ ] H J ) , GAM ( "t> \J ) 9 Ga-v„\ ( x -yj i 
1 A C l , X3 1 7 ) , 7 ACH { 7 > 2 f-C ) 

, T F A M 

{ o ) , mop ( -a ) , Dr LXF ( 7 ) > ANC- I 7 ! » ANGN ( ? ) 
I p D , > pu , j 7 v op , p AM » J FLO > ■ U~> > A L ? " * ~ 
R v X L ” , E P i H ; E A - 5 ? r P G » A A ! > 
vUP ( i 1 ? DU° i *J 


X ( -7 , l - r . ) 

X ( ^ ; >? 5 ( - ) » 0 5 N ( ^ 
5 v ( -< ) , hsn ( •> i 

* GAMA 5N ( 7 i 

5 0 ) * X i 2 3 0* >X N ! ] 5 0 
r o ) * t K ( ] A i *- / 1 

f", A 

1 

• W ( I 5 0 
» EM ( I 5 : 


1 


• \ ■ f . r A x 

:\ > i r- ^ i 


•‘.N (1=0 


» X J 


• < I \ F » I M r > N 


* JCH 


» <\ 


^ i\ H O U P ' r ) 


) = 


('■>)•> OUH ( 7 ) y i 

T A N ( QX^0X7«QXA ) 

QZ 1 * 7 AN ( GZ 7 ) +0Z 2 *Ta.Y ( OZa 


H ^ ^ ! 3 ) ' A L A* w* A ( /’’ » 7 ) 
) , XUP ( o ) » rvijp ( 7 ) , V.-PlJP ( '< ) 


T I? n F. = T >7 

T = 

STOP- A o W ;V S T P F A M doac 

M 'j - ? o 

Y T N 1 T - v M ( m 2 1 
v n c _ y ( a; 2 1 
X 7 .S = X ( N 7 ) 


: q T I c S 


o 


X I M 7 = XDS^O-LX- ( I R ] 
CO TO ( ]. *2»1 ) » TlYPc - 


i 

X-M'JSN ( t ) =XMUUP ( I PU i 

/7y 


oSN ( I ) =oijP ( T PU ) 

MM 


then ( 7 ) =pup ( t pu] 

„ 


OS M ( I ] =OUP ( T RU ) 

o 


T s M ( I ; =TUP( TRU ) 

C M5 A; ( 7 ] = ~MUP ( T PU ) 

o 


G A v. A S M ( T ) =G. A M A ! JP ( I » U ) 
Xr/WS.N ( 7 ) = X M U P ( I RU ] 
Pmos;,! ( T ] = P M 0 U ° ( I RU ] 

o 

10 

mSM ! I ) =^U D ( T PU ) 

AO ]A J = ] , M F D 

alpsm ( j , t ) = alpup ( j , i pu ) 

o 


y p = y p A. c 
p c T i ) R *• 

2 

GO TO ( pc , AC ) , T C L0 


•ap 

T p = IPU 



on TO A 1 


A A 

T R; 1 


A i 

M 7 = U 0 D (TP) 

p 


N o l = N L O ( t P \ 

r)C|_X = pcLXc ( TP ) 



GO TO ( S , FA ) , y rLO 

D 

73 

CO TO ( => , 6 ) , T UP 

A 

GALL cl t p ( 0 , NPL » N? , X M S V ) 

O 


call STcj ( o r LXc ( I P ’ , mpi_ ’ N7 ) 
7 c: ( jr i-ja-v, # trQ , 0 ) r, Q to fO 
a.O =■' ] L = N ° L , M 7 

O 


X m ( !_ ) = X ( L ) 4-priXc ( rR) 

Yp ( L ) =Y ( L ) + T ANj ( Th ( L > < * 0 E L X F 
aX =SOR T ( ( X,\ ( L ) -X ( L ) )«nl YW 


! YN(L)-V ( l_ ) )*#? ) 
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D 


j 


-.j 


O 

o 

o 

u 

O 

o 

o 

o 

o 

o 

o 

o 

o 


o 


o 

o 

o 


»X'jp( i a ) s x i 


6 

77 a q 

33 
26 

32 

34 
2 ? 

7 Q 0 


) GO TO ? y 


?0 

7 

7 9 


T W Pv = 1 

GALL rwCVO ( T VvP , r^x > !_ ) 

C O w 7 T N ! i c 
-0 TO c -2 

no s 6 |_ = N D L * •'-! 2 • 

p,n ~ =■ j= 1 , \icd 
W rvoT ( J i L 1 =0 . 
w^or\ ( j , l ) = 0 . . 
f'A r*-’ ( J . L ) =o . 
rfTMT } ,v,i;r 

yoi„n = y ( «7 ) 

Xijd ( jp] = x { N 7 ) + o E L X - ( IP) 

W R J T £ ( a » ■? 7 p o ) ? R , J <7 U • N 7 » N 

r a L L -LOW 

r a i_ L w - f. u 

T ^ ( A U I* i t 7 j — A T t } ^ A 

Y I m I o = ( X U ° ( JR ) — X D 3 ) # T A Nj ( FA v. ■>;- 4 \j G ! I X 0 ) ) — 

W R ! T F ( 6 > 7 7 4 o ) I R i I Ru » N ! 7 » N P L » XU P I I R ) , F A 7; 
F O R m A i 14,^1 ) 

YOLO = ( XOLD-XOS ) *TAN ( F Av*A : \6 ( I RD ) ) 
r O:'|T T M'jc 

! - ( I F L O , c O . ? ) r -n To 73 
rn To ( 28 »?2 » ? f 1 , t oij 
00 ? 6 JZ = NPL » n !2 
IF ( YOLO . L - . Y(JZ) 
rn\jT t vi jc 
FALL c 7ROI>(U) 
r-0 TO 77 
AO 74 = MP|_ <• M 7 

i- ( yolo ,g~. y(jZ) 
row! T NI.JF 
GALL FPROPilF) 

W = JZ 

K .8= J U - 1 

WR ! T p. ( 6 ,7790 ! Y I IX T 8 , YOL O 
copva.7 ( 8( 7XF14.6 ) ,718) 

Y L = 1 

■ R A T = ( YOLO-V ( <S ) ) / ( Y( IU) -Y( ■<;.• ) ) 

D5.N (<L5=P (KS)+R/t.T*(? (jU)- 0 C<S) 
qS,Y ( <L ) =0 ( <S ) + R,\T* ( 0 ( IL5-0 US) 

TS.V(nL) =T (.<S)+RaT*(T (jU)-T US) ) 
RHOS.N ( <L ) =Rh'0 IKS) J-RAl* ( RHO ( TU5-RHO ( 
THSN(<L) = T h ( '<s I +OAT “- i Tw ( I U 5 — T ;s (:<S) 
FY'S\ ( <L ) = F |V ( :<5 ) +PA.T* ( f.v { I U ) — ~ M (<5) 

G A v S 'V ( •' 1 . ) - G A v ( l 8 ) r.> A T *- ( G A v ( I 'J ) - 6 A M 
X N* W S M i < L ) = V.‘ (<S)+ X- A T * ( W ( TO)- '/• 

H.S.V { K l ) = HX (KS)h-RaT- (HX (IU)-HX { <S ) 
nO 13 j-] , M 8 o 

A L P G \ ! ( J , K L 1 = A L P ( J,KS)+SAT*(ALP < J -> T U ) 
row T 1 w: jp 

T F ( IF LO .- 0.7 ) r ,0 To 70 
r -0 TO IQ), TRV 

rO 7 JZ=npl,M2 

IF ( Y J MTS ,L". YN ( JZ ) ) GO T 0 7 1 • 

rOMT'i \*f if 

GALL fqqop(12) 

r-0 TO ' J ] 

r>f) 17 JZ--'-!PL,v? 


V.oc 
. Lv C) 


+YD8 


) GO TQ 77 


Y (JZ) »JZ».<S»tJ 


) 


<£) 5 


) 

»' < S ) 
US) 
) 


-ALP 


V_.y 


J«.<S) ) 
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o 


o 

o 


o 

to 

o 

o 

o 

o 


o 


.0 

-o 

o 

o 


If-' (VTMTS ,gp. Y n ( j Z ) ) GO in ^ j 

1? row 7 T NIK' 

r all cop or ( 1 3 \ 

31 T ! I = J Z 
X c . = t i.J — I 

WRIT :(?,i07CiC! Y I N T 5. • Y J ft 7 , Y G L P » Yft ( JZ i » Y 03 ? XOLQ » JZ * XS » I U 
0*7 <?0 "oov.aT ( M 7 X F 1 A . f. ) »?.!=,) 

R a I = (YtivTS-YN(nS) ) / ( Yftf US ) ) 

K L = ? 

dSN ( <L ) = D R ( <S ) -j-RaT* ( Pf\; ( j J ) -PR ( xS ) ) 

0 5 w ( < l ) = 0 ( •< S ) + P f T » ( 0 N ( T U ) - Q ' ( x S ) ) 

T S -v ( ,< L ) = T ft ( :< S I 4 . R a T * ( T M ( T U ) - T N ( < S ) ) 

dhOSN ( XL ) =RHON ( X S > j-RA T ( RuQN i Tu 5 -9H0N i <5) ) 

7h$mC<l) =T«N ( <5 ) +PAT* ( THN ( I U i -7-»\' ( KS ) ) 

FMSi\(<L ) = FMN ( <5 ) 4 -paT* ( "MM ( ? U ) — F M ft ! X S j ) 

Q A M A 5 ft ( .< L ) = G A Yft US I -+-k AT* ( G AA'-ft { i L : i - G A M ft ( <£ ) ) 

X Viv'SiV ( <l )' = Wi\m',<S) 4 -RaT*( '.v ft ( I G 1 — v.iY C <5 > ) 
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Dmudn = DTDN*V:SC*( 1.5/"/V( I 1-1. / ( T M ( I 5+19S.6) ) 

V I SC = ? « 2 7 F - o * ( 1 . 9 -* 7 M ( I ) ) **1 . v i ) . S*Ts\i ( I ) + i 9& . 5 ) 

VI SC = V! SC^-r ACTOR 

C 0 M O - V ' SC -“- C p X m ( I ) / PR 

n!JM=.0 

no 16 J=1*N S D 

16 DUV = DUv + AL D N ( J» I )*0C 2 ! (J) 

rv<n,\j= oMUD'N*CPXN ( I ) /°R+V T S C*DT;0\ ! / D R*DUM 
9090=0, 

SUM 1 =0. 

6 U m 2 = . 0 
S!JV0= o 0 

7 PRM 1 =ohOM ( T ) -“-CM ( I i 
AC 2 = X L r / D R-" V J SC 

IF ( ASS( YN ( I ) ) .LE. .GOOCl) GO 70 13 

AC1 = XL' r /PR*COS ( Th.M ( I ) ! *V I SC / Y.\ ( I ) *X J+XLE /PR*D.v.UDN 

C-.O TO 1 c < 

18 no 2? J=1»N-SP 

DADM >' J ) = ( ALPN ( J » I + i 5 -"- R A 7 !0l -ALPN v J » I ) * ( R A 7 I 0 1 -R A 7 1 02 ) -AL PN ( J » l - 1 ) « 
1 *o a 7 T 07 ) /SUV 

07 AON ( J i = ?.-!.-( ALPM ( J» 1 + 1 ) -*DELYi /SUM- A L PM ( J * I ) +AL?:\ ) ^0 EL 

1 Y? /SUV ) / 0 F L Y I / D :: L Y 7 

SO I N C J » I ) =TF,Rvi -!•- ( ALPM < J » T ) - ALP ( J, ! ! ) /OF LX- ( I R ) 

IF ( XJ.FQ. .0 5 S3 I N (J ♦ I J=AC?*‘D2ADN( J) 

OP PO- HI ( J I *WD0Tn ( J , I ) +H 1 ( J ) * S" T M ( J , I ) +OPRO 
S U r/ 1 = H 1 ( J ) * D A D N ( J ) + S U 1 
SUM? = h 1 ( J ) -“-D? ADM ( J ) + SUM? 

23 SUM? = D ADM ( J ) C D 1 ( J ) + S U m ? 

DD ! F= ( D MUON * S U M 1 D T DN * V I SC*SUM3+V I S C S U 2 ) / SCM 
T F R M = RhO M ( I )*CN( ! )->CPXN( I )*TN( I ) 
r, 0 TO 2 A 

19 no 11 J-1»NSP 

O AON ( J ) = { a L P N ( J » ! + 1 ) -^R A T I 0 ! - A L P N ( J » I ) ‘ : ' r 

1 ( R A T I 0 1 - R A T I C 2) - A L PN ( J ’ I - 1 ) «-RAT T 0? ) /SUv 
07 A DN ( J ) = 7 . * ( AL-pm ( J » I + 1 ) "-DF LY l /SUV.-ALP.N ( J » I ) + ALPN ( J» I -1 ) -DEL 


53 



o 
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o 


o 

o 

o 

o 


o 

o 

* 

o 

o 

o 


w 


o 


o 

.o 

-O 

o 




C 2 


. • » «-o,\ i : )■ 


1 Y ? / S U M ) / n " L v 1 /ri^LY? 

SO I N ( J » I ) = AC 1 *DA0N' ( J ) -r4C2*D?ADN (J) 

dpro= hi « J) ”WdoT;V( j, r )+hi ( j;-» j»i h-cpro 

SUM1=H1 ( J)*DADN< J) +SUM1 
SUM 2 = ( J 5 +SIJ : '’? 

11 SUM - = C-a DM ( J ) irr=>? ( J ) + SU'-' - 

DO I F = ( 0 v U D M * S U M 1 - D T D N * V ; S Q « 5 U V + V t S C*-SU M ? ) / S C V- 
TPRY = RhO\ ( T ) -::-eU ( J )sf PX.M( T ;-*TN( ID 

D?M ( I ) = ( S ! M f T H v «’!))/( Y* i Uifc a ! ~DUDN*VI SC+A . *CGS » T KM ( J ) ) -CM ' I ' -V 
1 *DTH,D\i/( YNi I )-»?. ) -V I SC-QM ( I 1*5 IN ! T HM ( I ) )*CGS(THN ( I ) ) / Y N ( I }*• 
2. *X Jt a k «V 1 SC^OUDN! vOT HDN / 3 * +/i . *V I DC-ON i I W?THDiVO 
1 # D T HDN /■».-?. *0N f ! )* S I\> ( 7«,v ( I ) ) a-DMUDN / ( 3.*Ym ( I ) )*XJ 
S 4 M ( 1 ) =DDIF+COND*COS(THN( ! ) ). *DTDN/ YN ( I ) * X J + CGN 
1 +4 „ -::-V I SC ( I ) **-.? /■* . * ( DT hD\*«*?-5 I N ( THN ( I ) J *D7HDN ) +V I SC*DUD:\-**? + 
? V I SC&QN (!)«•( COS ( Th-\J ( I ) ) *DUDN/YM ( I ) *XJ+D2U0N 5 +0N ( I ) -DUD.Y - QYUDM 
DIN ( I ) = V I S C * C 0 S ! T h N ( I ) )*DUQN/Y.N< I ) X J + V I S C * D2U0N+DUDN * D M U D N 
FDI SM ( I ) =S4N ( I ) /TEPM-Dl N v I ) * ( CPXN ( I )*TN ( I ) +QN ( I ) **? ) / ( TERM-OM i ; 
1 -D D RO/TpRM-S I M ( T HM ( r ! ) /YN ( I ) *X J 
C WR I 7E ( 6 » 1 CO.? ) D1M( I ) »S4M( I ) »DPRO»FD!SM( I ) » TERM > CPXN < I ! * I 

r-O TO 12 
24 roMT I MU C 
D 2 N ( I 1 = .0 

DlMU) = Vl$r*D2Un\! 

c D1 N ( !) sTERM} * f ON ( 7 >-0 ( I ) > /DFLXF { JR ) + ( PN < I ) -P { I ) ) /QELXFf I R ) 

S4\>{!) = D D I F + C O N D * D 2 T D N + D K D N DTD M 

1 +4.*V I SC *0N ( I ) * ( DTHDN«-S 2 -S I N ( THN ( I ) > *DTHDN ) +V I SC*DUDN**2+ 

2 V I Sr* ON l 1 )-> D-2UDN 

C S4N ( I ) = IFRM1 * ( hT.N ( I ) - Hi ( J ) ) / DFLXF ( IR ) 

FDIS:\! ( n =S4N ( I > /TEPM-DI N ( I I * ( CPX N { I )*7,\( j ) +QN ( I )**? ) / ( TERM*GN ( 

1 -DPRO/T PRM- . 5* ( ON ( II -Q ( I ) ) / ( TAN ( XMUN (in **2#QN ( I ) *DELXF ( I R.) ) * 

12 CONTINUE 
R c TURN 


; r 
) 


: - > 
X J 


2 CONTINUE 

DO 3052 IR=1»3 

I F ( I P .GT. 1 .AND* I TYPE .EC. 3> RETURN 

M01=ML0( JP 1-1 

NO 2 = MLO ( IR1 + 1 

N°1. = M0D ( r R ) +1 

NP 2 = N0n ( J R ) — 1 

T ( MQ 1 ) = T (NO?) 

O ( NO 1 ) =0 ( NO 2 1 
T m ( NOT ) =TH ( NO 2 1 
N4 = NLO ( r R ) 

M G = N 0 P ( T P ) 

Y ( NO 1 ) = ? . * Y ( N A ) -Y ( NO? ) 

Y ( NP 1 ! = ? . * Y ! NR ) -Y ( MP ? ) 

PC 5 0 J = 1 ? N S p 

A L P ( J » n O 1 ) = a L D ( J > N P 2 1 
5 0 A L D (J»M011=AL d (J»N02) 

Q ( NIP] ) =0 ( WP2 ) 

T ( vdi ) = T ( NP? ) 

Tw ( MP1 ) =Tm( no 2 1 

DUV=0 „ =- ( PWO (Ni)*0 ( M A ) +RHO ( NP ) *0 (NS ) ) 
DO 31 IT=NA»vp- 


= t T 

JM= I T-l 

IF (OHO ( NA ) so ( M A ) . GT . DUN ) v = NjR 


3 -A 


) 



J 
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o 
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o 
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o 
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o 

o 

o 

o 

o 

o 

o 

o 

o 

o 
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o 

o 




I F ( RHO (NA',+0 (\IA) .GT.DUmJ lV = NA +1 
I F { RHO (V)*0 )-DUM) 6 1 » 5 ? » 5 ? 

51 rON'T T NU r , ' 

^2 rhalf=y (?/)-(y (v,)_y (iiV))*(rhc (vn*o (M)-oUM) / (R hO (M)*q (y,)-RHO 
i ( i v ) ( im)) -y(na) 

RHALF=A Q S { PHALLI 

Vise =FACTOR{ IR)-»PHA.LF*PH0<\ ! A)*0(NA) 

Vi sex ( T R ) =v T SC 

W R T T E ( 6 » 1 0 0 1 ) VISC 

n'< ON = .0 

R^ll^N=.0 

no 13 t=\'a,iVp. 

npi_Y i =Y ( i ) -Y ( I -1 ) 

nFL Y2 = Y ( T + l ) -Y { T > 

PATIC1=0FLY1 /0Fi_Y2 
RAT T02 = nPLY2/0FLYl 

D'JDN = ( G ( I + i)*RaT I 0 i -0 ( I ) -> ( R A T 1 0 1 - R A T 1 0 2 ) -0 ( I ) *RA.T 10? i /SUM 
D2UDN =?.*(Q ( 1 + 1 ) *DELYl /SUM-0 ( I ) +0 ( I - 1 ) *05L Yp /5U/ V i ) /QE L Y 1 /DEL Y? 
Dl ON =(T < I+l )*RAT ICi-T ( I ) -* ( RA.I ICi-RAT I 02 ) -T ( I - v } *R A T I 0 2 ) / SUM 
D 2 T D N =?.*(T ( I +1 )*DFLYi /SUM-T ( I )+f ( I -1 1 *DE L Y? / SUM ) /D E L Y i / Or L Y2 
DTHDN =(TH( I + i ) * R a T 1 0 1 - T H ( I ) * ( R AT I 0] -RAT 102 )-TH( I~1 )*RA.T'!0?) /SUN 
D2THDN =2.*<TH ( I +]') *DELY1 /SU.v-TK < I ) + TH ( I - 1 ) *DEl_Y2 /SUM ) /0~L Y 1 / 
lon_ v 2 

CALL T hFRM ( I » I 5 
cond=vt sc*c d x ( n /pq 
n.lJY=,0 

no 15 J=1,NSP 

15 OUM=DUY+AL D ( J* I )*0CP1 ( J) 
nnco=0. 

SUM] = .0 
SUM? = . o 

suy?=.o 

AC?=XL C /PR*V! SC 

I F ( ABS { Y ( I ) ) .LF. .0000]) 00 TO 21 

AC! =XLF/PR*COS ( T.h (I) )*Vl Sc/Y ( I ) *X J+XLE/PR*DMUoN 

r-0 TO 2 2 

21 26 J = 1 * N S P 

DAON(J) = ( A L D (J» 1 + 1 1*RATI01-AlP < J ’ I ) * 

1 ( RAT 10] -RATJ02 )-ALP ( J, I - 1 ) *RAT I 0? ! /SUM 
D 2 A 0 ,v ( J ) = 2 ♦ * 1 A L D f J > I + 1 ) *DEL Y ] /SUM- AL D (J>I)+ALP <J*I-li*DEL 

1Y2/SUM) /OFLYi /OELY? 

S 3 t ( J » I ) = 0.0 

I- (XJ.FQ. .0) S 3 I ( J.» r )=AC2*D?ADN(JJ 

D D R0 = HO ( J) *',\’DOT ( J , I) +Hi<J5* S °! ( J ■> I ) + D?RO 

SUMi =K] < J ) #Dfi ON ( J ) +SUM1 

SUM?= HI ( J ) *[>2ADn (J)+SUM? 

26 SUma=. o A n m ( J ) -k- r o 1 ( J 5 + S U 2 

DD I F = ( D.YUDN*SUM 1+DTRN*V I SC^SUMA+V I SC*S*JK2 ) / SCM 

ter m. = r h o ( i ) *o ( i ) •> c ° x <n*T in 

r.Q TO 2 5 

22 no 10 J = 1,KISP 

OAr)N(J) = ( A L p ( J» . 1+] )*RATI0?-ALP f J> I)* 

1 ( RAT 101 -RAT 102 ) -ALP U> i-i )*RAT]02) /SUM 
DrADv(j) = ?•■»( A L P ( J » I + j ) -»DELYi /S'J.m-ALP (J»I)+AlP (J»I-1>*DE L 

1 Y? /SUM ) /p r LY 1 / 0 E L Y 9 
S3 I ( J » I ) =AC1*DADN ( J ) +AC2*D2ADN ( J ) 

DP, R0= H] * J ) +WDOT ( J , J ) +Hl (J)* S3 I ( J » I ) . +DPRO 


n 
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SUM 1 = H1 ( J ) :: OAON ( J ) + SUM 1 
SUM? = H ] ( J ! * 0 2 A 0 N ( J 1 +SU'-'2 
10 ,SUMM = . D*DN < J)*-CPi ( J)~SUM3 

DO I F= ( DMUD.N#SUyj +DTDIV-»VI5C*SUM OJ -VI.SC*SU.M2 J /SCM 
lrpw = RHO ( I ) -so ( I )*C°X (‘7)*T (I) 

07 ( I ) = ( S J N ( T H C I ) ) / C Y ( y ) * . ) *OUON* V I SC + 4 . *COS ( T’H ( I >) *G t I ) 

1 *DTHON/(Y ( I >*3. )-4.*VISC*-0 (T)*SIN(lH ( I 1 )*COS(7K ( I ) ! /Y (I 

2 ) *X J + 4 . *VI SC^DUDm*DTHDN/ A SC*C ( I) 02 ThD.N / 3 . + A • *Q ( I 1 *D 

? *DTHDN/3 .-? o*Q ( n* Sli\'{ Th i I i (3.*y ( I ) )*XJ 

S4( I ) =DOIF+COWD*COS(TH { I > 5 *D?DN/ Y < I) *X J+CONO*D2 TDN + DKDN* 
1 a.*VISC*G ( I ) **2 / ’. . * ( DTHDN** ?-S ! N ( TH ( I ) ) *DTHDN > +V I SODUDN** 
? V T SC*Q < I ) * ! COS ( Tm ( j ) ) *nUf>N/Y ( I ) -*XU+D2UDN > +G ( I ) «DUDN*CMUD 
OKI) = V I S C * C C S ( Th ill )*DUDN/Y ( I ) *X J+V I SC02UDN0UDN0KUDN 
PQ I S ( ! ) = S4 ( I ) /TERM-D1 ( i )* (C = X ( I ) * T (I)+0 { I )**?}/( TrRM*0' i 
1 -O d RO/TfRM-SIN(Th ()))/Y (n-x-Xj 
r - 0 TO 13 
25 CONTINUE 

TcRMl =rhO ( I ) *Q ( T ) 
o? ( n= .o 


) -* 
MU 


/ ~r 
i\ 


01 (M = VlSr*02UON 

S4 (I) =DDIF + C O N D *• D ? T 0 N D .< > .v *• D T D \< 

1 +4 . #V I SC*Q ( I ) *-*2 /"*•#( D T H DN np—S I N ( 7 H ( I ) ) -> DTHDN J +V I SC*DUDN**\~- 

2 VlSr*0 £ T ) *• 02Un,\> 

FDIS ( I ) =S4 ( I ) /TFRW-Ol ( t ) * ( c D X ( J > T (I5+Q (])**;?}/( TERM:*Q (J) 

1 -o 0 ro/t=-rm 

13 CONTj NU" 

30 52 CON T J MU~ 

PC TURN 


7 001 FORMAT (5H CAL2 * 5 ( 2XF1 h .6 ) >15) 

1002 FORMAT (TH COMPS 2 » 6(2X514.6) *15) 

1003 FORMAT (7(2Xci4.6) ) 

1004 FORMAT (BH COMPS 4 » 7 ( ? XE 1 4 . 6 ) » I 5 ) ' 

1005 FORMAT ( 5 ( 2 X F 1 4 . 6 ) ) 

1007 FORMAT (7h COMPS 7 ,8<2XF1.4.6) 5 

1.008 FORMAT ( 7H C0MPS8 ,4(?XF14.6> * 

C N0 



c 


o 

o 

o 

o 

» 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

p 

-O 

o 


• V - 

COMMON vr (200 1 »V» ( 2 CO) *XST.A(<= ) *NPT’(f») »I I < 6000 ') 

COMMON P ( 2 00 ] ♦ T ( p CO ) »L*( ?00 ) , T H ( ?0C ) , fM, ( ? 00 ) , Y { ?0G »Y? ( 200 » 3 ) 

COMMON ALP ( 7 , ?CO ) .Yi ( ?C0 ,5 ) >Y? ( 200 , 5 ) » Y3 ( 700,5) 1 YM ?GC > 5 ) 

. COMMON Y5 ( 200 , R ) , Y6 ( 2&0 ,5 ) ,Y? ( 200 , 5 i ,Y3 ( ?0C ,5 ) , Y9!2CC>5) 

common* fr AS* mcUPV 
EQUIVALENCE (OY,5CL v ) 

POUTVAL c NCF(OX,SCLX) 
lOO FORMAT ( 7 c 1 0 . 6 ) 

I 01 • format ( 6 F 10 .L) 

!05 FoqVAT(ftM 

c i c a S s 1 SP-r IF plot 

r I C A S = 2 FLOW P R 0 p - R T T fS PLOT 

C • NCJRVr NO. OF CURVFS PFR PLOT (X STATIONS) 

CALL PLOTS ( T T * ftCOO 1 
CALL PLOT ( .0 , 0. ‘ 

T conttnij- 

RFAD(S,iOO) ACAS,ACURV 

T C A S = A C A 5 

NCURV=acURV 

O0 22 LC=],NCURV 

R L * AO i , I CO ) XST A ( LC ) vAPT 

N P T ( L C ) = A P T 

r moT = NO* C~ POINTS IN- PROF! L c 

N P = N P T ( l_ C ) ' 

WR I T E ( 6 » 1 1 0 ) X S T A ( l C ’ 

II 0 FORMAT ( ! .0 X » 2 3 H PROFILES AT STATION X=,F 10.4) 

OO 12 T=1,NP 

RFAO ( 5 , id ) P ( I ) » T ( I ) *U ( I ) ,TH( I ) , f,m( I ) , Y ( 1 ) 

WRITE (6 *3 on P E I ) » T ( I ) » U ( J ) > T H ( I ) » EM ( I ) > Y ( I ) 

12 FONT I NU r 

PO 21 J = 1 * N ° 

RE AD E 5 * ] i 0 ) ( ALP ( K » J ) , i< = 1 ,7 ) 

C W R I T E ( 6 » 1 0 0 ) ( A L P ■( K , J ) » K = 1 » 7 ) 

21 CONTlNUF 

oQ 2 2 J= 1 « N D 
Yl( J»I f)=aLP{A»J) 

f vi=m2 

Y 2 ( J » L f ) =aLP ( R , J ) 
r Y 2 = 0 2 

Y3 ( J »Lf ) =aLP (3 » J) 
r Y 7 = H20 

Y4( J , Lf ) =aLP ( 7 , J } 
r Ya= N 7 

Y5{ J*LC)=P( J) 
c Y5 = D R- SSURf - P 

YS(J»LF)=T(j) 
r YA=TFMPFRaTORf _ T 

Y7( J,LF)=U(J) 

C Y7=VFL0CTTY - U 

YR ( J , L C ) = T h E J ) 

C Y8= FLOW a N c L F - Tm 

Y 9 ( J » L o ) =FM { J ) 

C Y9 = M A f m NOMA r P - cm 

Y D ( J , LC 1 =Y ( J 5 

22 FONT I NO f 


/* 


0 

r 

O R 0 F R A M r - P ( 

f u C p 


F A L L p A p T 2 




ro TO l 


0 


C N| r\ 




o 

o 

o 

* 

O' 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

.o 

o 

o 


SUpROOT TNf DAQi? 


common 

COMMON 

COMMON 

COMMON 

COMMON 


t r , 


VV ( 200 1 i YM ?00) »XSTA.(5 ) >i\PT < 5 ) » I I (&OCC) 

P { ?0C ) , T (?00 ) ,U { ?00 ) »Th< 2 0 0 i ,c;/,(jCCi >Y( ?.C0 ) ,YP ( 
ALP ( ■> >2 00) » Yi I ?00 »«0 »V;M ?0C» 3 > »Y3 ( ?00 , 5 ) * Y4{ 20 
Y 5 ( ?00 > 5 1 » Y6-C 200 » 5 ) * Y? C 20 0 r 5 ) > YP, ( 200 , = ) » Y9! 20 
»S» NCURV 


' a ; 
C » 5 ) 


100 
1 01 
I 05 


111 


2 5 


24 


FQU T VA.LF NCF ( DY » 5 C L Y ) 

p oo: valfncp ( ox>sclx i 

04 T a O0NF/4W00NF/ 
c 0 P M A T (7F1C.X) 

FORM-'.T ( 6 F 1 C . A 
FORMAT ( 44) 

RF AO t 5, 100} X M » X 0 > .5 ! 4 X > Y 0 » Y M > £ I Z Y 

SCL X = ( XM-XC ) /S I Z X 

S C L Y = ( YM-YO) /SIZY 

WR ! ” F < 6 » ) 0 1 } OY » SCLY ? DX »SCLX 

WR J T 5 ( 6 > 1 02 ) X 0 > X l v i » S I 2 X > Y 0 > Y M » S I 2 Y 

Y] =H2 

CALL A X I S ( ,0».0>3H Y »-f,StZY, # 
CALL AXTS( . 0»«0»14H H2 -RACTION 
oO 2 4 L C = 1 > N C OR V 
WRITE (6 >111) XSTAfLC) 
format (iox>?3h profiles at sta 
1 T T c 0 / / ) 

NP = NPT ( LC ) 
oO 25 J=1 *ND 
Y A ( J 5 = Y J ( J,I_C ) 

YY ( J) =YD i J,LC) 

YA(NP+1)-X0 

YA ( \‘P~2 ) =SCLX 

YY ( NP+1 ) =Y0 

YY { NP-r2 ) =SCL Y 

CALL LINF(YY>YA>N D ,1>0,3) 

fONT 1 N U c 

X A 1 U = S r 2 Y j- 2 . 

YN1J=,0 

CALL PLOT ( Xi\U » YMU » -3 ) 

Y2= 02 

RF AD ( 5 > 1 00 } XM » XO > SJ ZX » YQ > YM > S I Z Y 

Sc L X = ( X M - X 0 ) / 5 I Z X 

•SCLY= (YM-YO) /SIZY 

WR I T F ( 6 » ]. 0 1 ) OY > SC.LY > OX > SCLX 

WRTTF(ft.lOi) XC>XM,SiZX>YC,Yy.,STZY 


> YO>DY) 

* 14 » S I ZX j 


9 0 • » X C > D X ) 


ION X= Fl0.4>23H ARE NOW BEIN' 


AX IS ( . C , .0 > 3H Y ,-3>S]ZY,.' 
AXI5(.0i«0»14h 02 FRACTION 
jpv 


27 


26 


CALL 
CALL 

00 26 LC = l»NCl 
m p = m p j ( l c 1 
oO 27 J=! »NP 
YA ( J ) = v? U.LD 
YV( J)=YD( J.LC) 

YA ( N P 1 ) =X0 

YA ( (VP-? } =$FLX 

YY ( NP-f X ) =Y0 

YY ( NP+2 ) =SCLY 

FALL LT,NF(YY »Ya.»NP»1 »0» 

FONT I NU-f 

XN!J = 5 T Z Y + 2 . 

VMIJr . 0 


, YO ,DY ) 

, I 4 > S I Z X > 


> X 0 > D X ) 


3 PLO 


58 



0. 

o 

o 

V 

o 

6 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

p 

-O 

o 

o 

o 
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TALL PL^'T ( , w .\iU , YNU , -3 ) 

R= AO ( 6,100) X i 1 ,.X0 »St 2X ,YC»*.Mr5 tZY 
SCLX= ( X y - X 0 ) / S I Z X 
Srt_Y= ( YVi-YO 1 / S I Z Y 
WR ] 7 E ( 6 ’ 1 0 i j DY , SCLV ,DX ,5CLX 
w R T i F ( A » 1 C 1 ) XO , Xiv , S i Z X , YO , YM , S I Z Y 
v ' = uj»n 

CALL A X I S ( • 0 , „ 0 » 5 H Y , - ^ , 5 I ZY , 

CALL AX : 3 ( .0 , .0, 1AW H 2 0 FRACTION » • 14 »SIZX » 90.»XO»OX) 
no 3SL0 = l«MrlJRV' 

M D =VPT (Lf) 
of) 7fi J=i,ND 
YA ( J>=Y? ( J,LC) 

YY ( J ) =Y p ( J,LC ) 

Y A ( N P + 1 ) = X 0 
YA ( NP-2 5 =SrLX 
YY ( N P + ] ) sYO 
YY { NP-2 ) = ScLY 

CALL L T N - ( Y Y » Y A » N P , 1 » 0 » 5 ) 

CON'TtNU- 
XNU = S • 7Y + 2. 

' v MU=.0 

CALL Plot ( XNU , YNU »-•* ) 

P-AO( 6,100) XM,X.0,SlZX»Y0*YM,Sr2Y 
Sr L X = ( XM-XO ) /Si ZX 
Scl Y= ( Y M - Y 0 ) / S I Z Y 
\u- N2 

CALL AXIS(.0,.0»?HY , - 3 , S I ZY , . 0 , YC , D v ) 

CALL A X I f ( . 0 » . 0 , I A +' M2 FRACTION , 1A,SIZX, 90. ,XO»OX ) 
00 36 L C = 1 » N C U R V . 

MP-MDj ( L r ) 
no ?Q J = i , N o 
YA ( J ) = Y A. (J.LC) 

YY ( J ) = V P ( J , LC ) 

YA C f’JP + 1 ) =X0 

YA ( NP + 2 ) =SCLX 

v Y ( X) P + 1 ) = Y 0 

YY ( NP + 2 ) =SCLY 

CALL L ! w p ( Y Y . Y A > N P » ] , 0 » 3 ) 

rOVTTNU- 

XNO=StZ v +2. 

VMM- f 0 

'•ALL PLOT ( XNU, vml;,_3 ) 

R ^ A 0 ( 5 , 1 0 0 ) XM , XO , S I ZX » YO ! YY, » S I Z Y 
SCLX= ( X M-XO ) /S I Z X 
5- <" L Y = ( YiA-YO ) / S I Z Y 

Y 5 = D R cr 5 c, \j p c _ d 

CALL A X j S ( . 0 , . 0 , 3 H Y , - 5 , S I Z Y , . 0 , Y 0 , 0 Y ) 


5 u 


CALL AX J S ( . 0 , . G , 1 AH 
00 3? LC =1 * NCURV 
NP = N c T ( LC) 
nO 3 0 J = 1 , N P 
YA ( J ) =Y=. ( J,LC ) 

YY ( J ) = Y P ( J, LC ) 

Y A ( M p + 1 ) = X 0 
YA ( NP + 2 ) = Sf L X 

Y Y ( N D + ] )=Y0 
YY ( NP + 2 ) =5CL Y 


■CC'H- 


1 A • S I Z X 


»xoox ) 
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o 
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o 

o 
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CALL L ! N F ( Y Y » Y A * N P ♦ 1 >0i?) 

3’7. FONT I M'.ic. 

X N U = S t Z Y o. 2 . 
v N'l=. 0 

FALL PL'ni iXwli?Y!vJ,-^! 

R c Af) ( 3 » i CG ) X N- » X 0 » S T 2 X » Y 0 » Y V. » 5 T Z Y 
Sri_X= ( Xy-Xo ) /S i Z X 
Sri_Y= ( Ym-YO ) /ST ZY 
r Y6 = 7 fmoc-qaTL'R- - 7 

CALL AX ! S( .C , .0 » 3H Y »-3 » 5 I 7.Y , . 0 , vq , ^Y ) 

CALL AX I S ( , 0 , , Q . 1 4H T VpFR A TORE * 14*SjZX, 9C, .XQ»DX ) 
no lf = i ,Nirupv 
• MP=NPT(Lr) 
r-0 3 1 j=i .,\!P 

Y A ( J \ = V A ( J »' L C ) 

31 YY t j \ = YP (J,Lr] 

Y A ( M D + 1 ) = X 0 
YA ( NP-i - ?. ) =SCLX 

■ YY ( NP+1 ) =Y0 . 

YY { NP + ? ) -SCL Y 

CALL L T M “ ( Y Y * Y' A » N P » 1 » 0 * ?. ) 

3 8 rON'T T N ! j sr 

XMU=S:7Y+2. 

Y\l!J=.0 

call plot ( X\'U» YNU»-3 ) 

' READ I 5 > 1 00 ) XM»XO»5.lZX»YO»Yv,,Sl7Y 
SflX = ( XM-XO ) / £ I Z X 
■ SCL Y = (YM-YO) /Sr ZY 
r Y 7 = V f l O r T 7 Y - U 

CALL AXJS1.0*.0*3H Y ♦ - ? » -S T Z v * . 0 > YO » D Y ) 

CALL AX I S ( , 0 » . 0 » 1 4H VELOCITY , 14>SIZX» 90.,XO>DX) 

no 39 Lr=l»NCUPV 
NP = N p 7 ( Lr ) 

nO 32 J = ] ♦IMP • 

YA( J)=Y?( J,LC) 

32 YY( J)=YO{ J,LC) 

YA ( \tP4-] ) =XO 

Y A ( N P + 2 ) = S C L X 
YY ( MP-1 )=YO 
YY ( ,\jP + ? ) =SCLY 

CALL LIN!f (YY»Y'a*NP»] »0*3) 

39 CONTINUE 

XNU=S T Z Y + 2. 

Y M U = . 0 

FALL PLOT ( XN'U » YNL! » -? 1 

R" AO ( 8 > 1 OC ) X M » X 0 » S I Z X » Y C » Y .v * . c . T Z Y . 

Sfl X = ( X.v-XO ) /s I Z X 
Sr l Y = ( Yv,-YO)' /Si ZY 
r Y ?, = FLOW ANFiL- - Th 

CALL ax I S ( . 0 , .0 » 3H Y ,- 3 ,S I ZY , . 0, YO , r>Y ) 

CALL AXIS! .0, ,0» 14HFLC 1 // ANGLE * 14»S!ZX» 90. ,XO*DX ) 

00 AO LF=1-NCUPV 

MP=NPT ( LF ) 

no 33 J=i,NP 

YA(J)=Y.R(J,LF) . 

37 YY ( J ) = Y D < J » LC ) 

Y A ( NPj -1 ) =xO 

Y A ( HP-*-? ) = SC L X 


n 
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o. 

1 


• . . ' c 

0 


YY (NP- 4 -! 1 = Y C 

Y Y ( M P + ? ) = 5. r L Y 

CALL LTNF(YY,Ya»NP.1 ,0»3) 

o 

40 

ro^T T NU C 

i 

i * 


X N 1 ) = 5 T ? Y -i- 2 . 

1 


Y.\ll=.0 

j o 


rALL plot ( XNUt YNLS-3 ) 

I 

i 

| **■ 


R~A0 ( 5* ] 00 J X Ni > X C » S I Z X > Y 0 > Y M » S I Z Y 
$rLX= ( X m - X 0 } /Si Z X 

I o 


5rLY= ( YM-YQ ) / S I Z Y 


r 

YQs M h H .V'JYRFP = eV 

CALL AX I S ( . 0 » . 0 » 3H Y , -3 , S ! Z Y . . C » YO » DY ) 

o 

i 


CALL AX 1 S ( . 0 , .0 , 14HMACH NUMBER .♦ 14*SIZX» 90.,XO»DX) 

^0 41 LO= 1 »ncu?v 
ND = \'PT ( LT ) 

o 


no 3^ J=1,NP 
YA(J1=Y9(J»LC) 

o 

34 

YY (J)=YD(j,Lr) 


Y a ( NP+ 1 ) = X 0 

Y A ( N P - ?) = S C L X 
YY ( fjjp+l 1 = YC 

o 


YY ( fvjP + 2 ) =SCLY 

TALL L!NF( YY»Ya»MP»1 »0»3 ) 


41 

rONT ! N»J C 

0 


XN!U= S T 7 v + 2 . 

YNU= • 0 

CALL p L 0 T ( X N U » Y N! (J » - 3 ) 

o 


or/inllUO^ PUNN 

o 


T F { OOA' C . -0. OUMNn GO TO 50 
pc- TURN 

50 

CONTINUE 

o 


GALL PLOT ( 2. » 2. *999) 
. pc-turn 


PNP 

kJ 

O 

i 



o 



! ° 



o 

4 



* o 



o 

• 

• 

o 
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TABLE 2 Input Cards 


Card No . 1 


SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

IJ 

0 for 2D, 1 for Asigyiranetric 

1-5 

15 

IOUT(l) 

printed output for region 1 every 
IOUT step 

6-10 

15 

IOUT (2) 

printed output for region 2 every 
IOUT step 

11-15 

15 

IOUT (3) 

printed output for region 3 every 
IOUT step 

16-20 

15 

NSAVE 

not used 

21-25 

15 

JCHEM 

0 frozen flow } l finite rate 

26-30 

15 

ISTART 

0 initial profile input according 
to format I 

31-35 

15 


1 initial profile input according 
to format II 




see pages 66-68 



Card No. 2 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

NMAX(l) 

Maximum number of points permitted in 
region 1. If this number is reached 
the program terminates and punches 
results at the last station 

1-5 

15 

NMAX(2) 

same as above except for region 2 

6-10 

15 

NMAX (3) 

same as above except for region 3 

11-15 

15 


NMAX (1) < NMAX (2) < NMAX (3) < 150 


Card No. 3 


SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

XPRO(I) 

X-stations at which calculated profile 

1-10 

7F10.4 


data is to be punches -to be used later 
for plotter- _ MAX number 14 (2 cards) 

11-20 


Card No, 4 


•% 


SYMBOL 

DESCRIPTION - 

COLS. 

FORMAT 

PR 

Prandtl number - constant (1.0) 

1-10 

7F10.4 

XLE 

Lewis number - (1.0) 

11-20 


FACTOR (1) 

Empirical coefficient (K) in 
viscosity model for region 1 

21-30 


FACTOR (2) 

ii u ii ii ii 2 

31-40 


FACTOR (3) 

ii n ii tt n ^ 

41-50 



Typical value of K = .02 



Card No. 5 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

CB 

Number of coordinates on lower wall 

1-10 

7F10.6 

COWL 

Number of coordinated on upper wall 

11-20 


Card No. 6 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

XRC (J ) 

X coordinates of lower wall (ft) 

1-10 

7F10.6 


C values required (MAX = 100) 

11-20 


Card No. 7 




SYMBOL 

DESCRIPTION 

COLS . 

FORMS? 

R'RC(l) 

First y coordinate of lower wall (ft) 

1-10 

F10. 6 
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Card No. 8 


SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

THC(J) 

Lower wall angle at above x-coordinate 
(degrees) 

1-10 

7F10.6 


CB values required 

11-20 

• 


Card No. 9 


* 

• 

• 


SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

XRC(J) 

x-coordinates of upper wall (ft) 

1-10 

7F10.6 


cowl values required 

11-20 


Card No. 10 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

RRC(cowl) 

First y coordinates of upper wall (ft) 

1-10 

F10. 6 

Card No. 11 




SYMBOL 

DESCRIPTION 

COLS. 

. FORMAT 

THC (J) 

Upper wall angle at above x-coordinate 
(degrees) 

1-10 

7F10.6 


cowl values required 

11-20 


Card No. 12 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

EPP 

Not used 

1-10 

7F10.4 

EPTH 

Error criteria for flow angle in the 
shock wave point calculation 

11-20 


Card No. 13 




SYMBOL 

DESCRIPTION 

COLS . 

FORMAT 

ANG(l) 

Magnitude of the initial shock wave 
angle in region 3 (degrees) 

1-10 • 

7F-10.4 

ANG(2) 

II II It II II II 

" " "2 (degrees) 

11-20 


DYASL 

Maximum distance, Ay, between the shock 
point and lost mesh point 

21-30 
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SYMBOL 

DESCRIPTION 

COLS . 

FORMAT 

DELY 

Rot used 

31-40 


XBP 

Initial x-coordinate (ft) 

41-50 


FLO 

1 Diverging Shocks 

2 converging shocks 

51-60 


XSTOP 

Maximum x-station calculation will 
proceed to (ft) 



Card No. 14 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

PUP (2 ) 

Upstream pressure for region (2 ) (psf) 

1-10 

7F10.4 

TUP (2 ) 

Upstream temperature for region (2) 
<°K) 

11-20 

- 

EMUP(2) 

Upstream Mach number for region 2 

21-30 


GAMAUP (2 ) 

Upstream Gamma for region 2 

31-40 


XMWU P (2 ) 

Upstream molecular weight for region 2 

41-50 


Card No. 15 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

PUP (3) 

Upstream pressure for region 3 (psf) 

1-10 

7F10.4 

TUP (3 ) 

Upstream temperature for region 3(°K) 

11-20 


EMUP (3) 

" " Mach number " " " 3 

21-30 


GAMAUP (3) 

Upstream Gamma " " " 3 

31-40 


XMWU P (3) 

Upstream Molecular Weight " 11 ^ 

41-50 


Card No. 16 




DU P (2 ) 

Upstream flow angle with sign for 
region (2) (degrees) 

1-10 

7F10.4 

DUP (3 ) 

it ii ii ii it ii 

" "(3) (degrees) 

11-20 


DXN 

Not used 

21-30 


DIN 

Not used 

31-40 
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SYMBOL ■ 

DESCRIPTION 

COLS. 

FORMAT 

DXIK 

Not used 

41-50 


COALT 

Characteristics coalescence criteria 

on /\y 51-60 

i.e. minimum distance between mesh points 


Card No, 17 ■ 




SYMBOL 

DESCRIPTION 

COLS . 

FORMAT 

ALFUP(J ,2) 

Species mass fraction upstream of 
region 2 

1-10 

11-20 

* 

7F10.4 

ALPUP(J ,3) 

n ti ii ii it ti n 

region 3 

1-10 

11-20 

7F10.4 

Card No. 18 


• 


SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

ITYPE 

1 uniform flow in region 1 

1-5 

15 


2 nonuniform flow in region 1 



Card No. 19 




SYMBOL 

DESCRIPTION 

COLS. 

FORMAT 

NL0(1) 

Lower index number for storage of data 
region 1 

in 

1-5 

15 

NOP (1) 

Upper " " " " " 

region 1 

6-10 


NL0(2) 

Lower " " " " " 

region 2 

11-15 


N0P(2) 

Upper " " 11 " " 

region 2 

16-20 


NLO (3) 

Lower " " " " " 

region 3 

21-25 


NOP (3) 

Upper " " " " " 

region 3 

26-30 


Initial Profile 




FORMAT I - If ISTART = 0 the initial profile is read in as follows 
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Card Mo. 20 


SYMBOL 


DESCRIPTION 

• 


COLS. 

Y(I) 


Y coordinatss of 
(ft) 

data line in 

region 1 

1-10 



(NOP (I) - NLO(I) 

+ 1)/, cards 

required 

11-20 

Card No. 

21 




# 

T 

SYMBOL 


DESCRIPTION 



COLS. 

P(I) 


Pressure at above y-coordinates (psf) 

1-10 

T(I) 


Temperature " 

fl 1! 

" (°K) 

11-20 

Q(X) 


Velocity " 

tl 11 

" (ft/sec) 

21-30 

TH (I) 


Flow angle " 

11 11 

(degrees ) 

31-40 



NOP(I) - NLO(t) +1 cards requiied 


Card No. 

22 






SYMBOL 



DESCRIPTION 





COLS. 

ALP(1, I) 



Mass fraction 

of 

H at above y- 

•coordinates 




in region 1 





1-10 

ALP (2 ,1) 



II >1 


II 


0 " " 

II 

11-20 

ALP (3,1) 



fl (1 


II 


h 2 o 

II 

21-30 

ALP (4, 1) 



11 11 


11 


H 2 

11 

31-40 

ALP (5, I) 



»f M 


11 


°2 

11 

41-50 

ALP(6,I) 



11 It 


II 


OH " 

11 

51-60 

ALP(7, I) 



II fl 


11 


N " 

2 

II 

61-70 

Repeat cards 

21 and 22 

for 

regions 2 and 3. 



Input region 

1 

CARDS ; 

20 

followed by Cards 

21 ’ 


II II 

II 

2 

II 

II 


II 

II 11 

It 


II If 

II 

3 

If 

II 


II 

II II 

If 



FORMAT 

7F.10.4 


FORMAT 

7F10.4 


FORMAT 

7F10A 
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Format II (Computer punches data according to this format) If Is tart 
input is 
Card No. 20 


SYMBOL 

DESCRIPTION 




COLS. 

P(I) 

Pressure at y 

in Region I 

(psf) 

1-10 

T (I) 

Temperature " 

11 

11 

(°K) 

11-20 

Q (I) 

Velocity " 

II 

11 

(ft/sec) 

21-30 

TH(I) 

Flow Angle " 

II 

11 

(degrees ) 

31-40 


BLANK 




41-50 

Y(I) 

Y-coordinate 

II 

It 

(ft) 

51-60 

There are ^ 

(NOP (I) - NL0(I)+1) cards 



1= 

1 





Card No. 21 






SYMBOL 

DESCRIPTION 




COLS. 

ALP(1,I) 

Mass fraction 

of H at Y (I) 

in region 1 

1-10 

ALP (2, 1) 

11 II It 

0 

11 

11 11 

11-20 

ALP(3, I) 

II II II 

h 2 ° 

11 

11 11 

21-30 

ALP (4 , 1) 

11 II 11 

H 2 

II 

11 11 

31-40 

ALP(5, I) 

It II It 

°2 

II 

II II 

41-50 

ALP(6,I) 

II II II 

OH 

II 

II II 

51-60 

ALP(7 ,1) 

II 11 II 

K 2 

II 

tl II 

61-70 

Same number 

of cards as card 

20 




Input region 

3 Cards 20 followed by 

Card 21 


THEN " 

1 11 U 

11 

11 

21 


THEN " 

2 11 ii 

II 

II 

21 



1, the 

FORMAT 

7FI0.6 


FORMAT 

7F10.6 
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Table 3 


Sample Input Cards 


■tr j 1 


. 3?S 

■' / .*i 

> »' • J 

. CcOul 


.ZHO. ! S 
X'!.G i! 


• J ■£■ 

. 10 

• “ - 

. 14 

a L. 

• o 



I ■ 

1 

— 11. 

— - « 

- 1 X 

* V* t 

-13. 

—1 0. 

■ 1 

-G • * 


1 

1 1 • 0 
1 

£, 0 

• J • u 

4 • 1 J 




! 

i -o... 


» u 

\ 0 




i 






. 


1 

1 . Ouuul 

» L ; U*J x 

« UL«J i 





1 

£ j .3 



, ?' ?’:■ 3 

c. a 

.5 



1 | ^5; 

Ot i -j 

1 • 0*vC 

£0. 04 S 




J 

-PSOi? 

o o 

* u. * •- * * ‘ * 

, , , .O' : 0 . 

* * * * > *_r*v*— • 

■ •'*•••’£? t?x-9 - 



_ 


| —o ■ 

;■ r- -J c- 

• OUiy.‘ 

62. • 

• COG GO 

0. 0001 J 




SI - 1 - 

• V 

• O 

« £_•£> £- 

% ecu 

• i* OO 


1 

1 

! i .o 

; i 

.0 

• C-Oil 

. uL‘U 

a 700 

■■ 


1 

t 

ii 


A •> 



. 

. 

1 


’ |2 1 1 6 . \ 
|2116. 
2116. 
2116. 


1050. p- 

|05O. 

•050. 


pi 3.-. 

|r- 4 o? 

• 0 
. 0 
.0 


21 16.0 
£1 16,0 
2116.0 
£1 1 6 .0 
2116.0 
j. 6 I 1 




'£• 200 . 
-j5£00 
'£• 200 . 
• |5200 
£ 200 . 


. 0 

:i-c* ' . 

|.o 

1-3.0 


11.5 
11-7 
ji . 0 
P-2 r . 

. Gob 
> •_> '-"J 
. COG 
110. 
-p,?, 
£.6 


...... 



Card 20 

5 0 

Table 3 Sample 

Input Cards Region 3 




o 

P 

• j. • > 1 <> ■' n i •- 

i ! ‘‘l 

T 

■ . 4 7?? r :'-iXi3 

! i-.'fdt r ; a t a 

u 

. 9625422*0 4 


Thc-ta 

. 3. 1 f) 6 1 6 2 + r, 2 

Y 

- . 1 4 1 li ’) 0 F - 0 I 

[ii 

• i 13 

: I ' t 

. a o p 6 5 1 r>:-3 

. 941072 F 4 . 0 4 

- 

. 3. 0 6 2 6 4 P. + o 2 

-.500 0 0 0 F - (, ? 

? 

• 1 59300! v 

! ’ *t 

. 6 40 84 5F + I;;< • 

. d 6 C vi 7 > 7 h + • 1 4 


. 3. 0 (j 3 a 7 <- + n d 

. 6200002-0? 

>r'.; 

. 1 6 2 2 3. (:■ ! - 

1: *} 

. C 3 1 I 1 'i r + i- 3 

. 752 4.1 3 OO 4 

- 

. 3 5 3 5 5 0 F + n 3 

. 2160002-0 3. 

? 't 

. 1 6 o 4 3 1 F + 

:> /. 

.161/ i.i r. + n 4 

. 6 6 3 0 9 6 ' .4-0 4 

- 

.317 '-oGF + fii 

.3450 0 OF-Oi 

5 s 

. 1 5 7 9 3 1 2 > 

». *\ 

. 2114692+04 

. 6 2 2- 4 0 72*04 

- 

.9767302+01 

.42600 G 2 - 0 1 

5*6 

• 1 6 0 '■ cl (.. i 


. ? 6 ICO o F + 0 4 

. 1 0 6 6 2 7 E +■ o 4 

- 

. 1 1 7 6 9 3 c + o 2 

.5000002-01 

57 

. 1037792 + 

0 '♦ 

.2 37 623 F -04 

. 564 4 4 9f + 0 4 

- 

. 1 2 o n 2 6 E + n 2 . 

. 581 0 OOF- 0 1 

5 a 

. 1 7 7 o 5 P + 

•j ** 

,294? i'i 1 [' * G 4 

. f> (, 4 3 o 2 L -> 0 4 



, 6 8 ?. ft 0 Q E “ 0 2 

5 0 

. l923''5!:v 

,* •? 

. 2 9 81162+0 4 

. 5 e 6 0 5 F -i- 0 4 

- 

.8927102+01 

, 7 5 6 ! i 0 0 F - 0 1 

6 0 

. ? 7 5 4 7 9 !; v 


. 30 17962 + 0 4 

. 5 3 6 0 6 2 F 4- 0 4 

- - 

. 7 4 7 7 4 0 c + o 1 

. 8 2 6 0 0 0 F. - G i 

6l 

. 2-695 ii + 

•: *\ 

. 2 I.i 4 5 4 7 F: - 4 

. 5266972+04 

- 

, 61p65 0 F -> o 1 

. 9 0 4 0 0 0 2- 0 1 

6? , 

. 2 7 6 7 1 1 / r + 

‘t 

. 4 0 6 ^ -■ 4 t; + 1 ; 4 

.5172372+04 ' 

- 

. 4903302+01 

. 9770 0 0 2-01 

63 

* 2 4 + o a |- r 

; . z. 

. 2 0 7 2 9 6 F. + 0 A 

-. 51 0<>l4F + 0 4 

- 

. 3 7 2 4 2 0 2 + n 1 

.10 t! ■-< i.i 0 2 + 0 C 

64 

. ? •- o 4 s i + 


. 2 0 6264 {->0 4 

. 5U5923F + 04 

- 

.2612602+01 

. 1 1 1 7 o 0 2 + 0 e 

65 

. 2 4 4 6 7 .6 !-. + 

■* 

.2094572-04 

. 6021562 + 04 

- 

.1592002+01 

. 11550 02 + 0 0 

00 

. 2 4 6 1 6 4 1: -r 

f- « 

. 2 1. 0 62 6 1'. + f; 4 

. 4 9 9 2 2 7 [" + 0 4 

- 

. 693-40 02 + 0 0 

. 1 2 p 1. 0 0 2 - 0 G 

67 

• 2 5 0 6 7 O F + 

»; ‘1 

.2121445+04 

.456 9 1 5 F + -j 4 


.63000 02- nl 

.1317 -J 0 2 + 0 0 

6 6 

. ? 5 2 6 3 6 f; + 

• * *i 

. 21327 0 F + 0 4 ' 

. 4 94 964 F+ 04 


.66 o 9002+0 0 

.13 5 2 6 0 2 + 0 f; 

6 0 

. 2 5 4 < i 7 7 F + 

r% ‘r 

. 3141 6 6 F. - 4 

. 4 93 2 37E + 04 


.1110902+01 

.1446002+0 0 

Vi) 

. ?9676 72 + 

L' ^ 

. 31434 ?R + 04 

. 4 V 1 o. 2 1 <-. + 04 


.142.1 4 0 E + o 1 

. 151 0 0 0 2+ 0 0- 

7 1 

. 2 5 3 5 5 0 F + 

fN 

.31529 1 F + 0 4 

. 4 3 9554 F + 0 4 


. 1 6 0 2 9 0 E + n 1 ' 

.15 /20 02 + 0 u: 

72 , 

. 2 6 0 (1 6 7 2 + 

• • *■'; 

. 31540 91-4-04 

. 4 6 7 6 4 0 [ : + 0 4 


.1699102+01 

. 1 6 5 4 0 0 F + 0 0 

73 

. 2 6.1. 1 69P + 

ft /t 

. 214oF2F-G4 

. 4 d 6 di 4 9 r - + 0 4 


.1732702+01 

. 16 9 4 0 0 E +■ 0 0 

74 

. 2 01 0642 + 


. 3 1 2 3 5 6 F - 0 4 

. 4 6 6 4 9 d F + 0 4 


. 1678202+01 

.17 5 4 0 0 2+0 D 

75 

. . 2- > 24;-6Fv 

ij *'i 

. 3 0 4 736 F + 04 

* 4 6 6 4 v > F + u 4 


. 1 4 3 6 4 0 2 + 0 1 

. 1813002+0 0 

/ 6 

. 2 6 2 7 2 5 F + 


. 2 3 6 3 9 2 F; -:• 0 4 

» 4 3 o 9 7 4 v. 0 4 


. 9 3 3 3 0 0 E + o !3 

. 1 8 6 7 0 0F + 0 0 

7“7 

. 2 o 3 0. o [■> 

. . * i 

. 2 64 0 71 p + i ; 4 

.490730 + v 4 


. 5 7 o 2 0 0 r + f, 0 

. 1 915 0 0 2 + C.C- 

7 6 

. 2 6 2 7 :v7Ft 


. 2 1 4 5" 3 F — 0 4 

. 4 9 5 9 0 5 F + 0 4' 


.1485002+00 

. 1 9 5 6 0 0 2 + 0 0 

7.9 

. 263274P v 

f; 'j 

. 3. 6 3 (> 6 9 F 4 0 4 

.5 0 2 4 72 F + 0 4 


.5935002+00 

. 2 0 0 1. C 0 F + 0 u 

3 0 

. 2 4 9 2 !.- 5 P. + 

n <• 

.lG235.6F.4-0 4 

. 5 n 9 Oco + ,i 4 


. 3. 6 1 3 4 0 E + n 1 

. 2 1 1 0 0 G E + 0 G 

8 3. 

. 2 34 55 01- -r 


. 9 313762.4.0 2 

. 5 1 4 0 39 F + Q4 


. 13p4 4 05 + nl 

. 2? 0 6 0 0 t: + 0 0 

62 

. 224 37 4 F-. 

0 'i 

. 9 1. 4 P 4 6;-. + 0 e 

. 6 j. 6 6 9' 6 r. — i) 4 


. 7774002+00 

. 2 3 0 3 0 0 F + 0 0 

i*i 3 

.2177212- 

*. # i 

.90695 IF 4 . 0 3 

. 5 1 c 3 9 P i-. + 0 4 


. 3 6 6 7 0 0 E + o G 

. 24 0 3, 0 G F + 0 0 

o * 

O *1 

.2-1417 6 F + 

r. 4 

.9029052+03 

.5193305+04 


. 1 4 6 3 0 0 2 + 0 o 

. 2 5 0 f: 0 0 E + 0 0 

35 

. 2+2526F+ 

0 4 

. 9 0 1 C 7 6 F. 4 - 0 3 

, 519 / P 2 ;:.+ il 4 


.4860302-01 

.260 G '-.i 0 r + OC 
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2b. DIVERGING SMOCKS 

Fig. 2 Simplified Combustor Flow Field & . No: 
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Figure 7 Characteristic mesh 



